Neutron diffraction study on phase transition and thermal expansion of SrFeAsF 
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The magnetic ordering and crystal structure of iron pnictide SrFeAsF was investigated by using neutron 
powder diffraction method. With decreasing temperature, the tetragonal to orthorhombic phase transition is 
found at 180 K, while the paramagnetic to antiferromagnetic phase transition set in at 133 K. Similar to the 
parent compound of other iron pnictide system, the striped Fe magnetism is confirmed in antiferromagnetic 
phase and the Fe moment of 0.58(6) p B aligned along long a axis. The thermal expansion of orthorhombic phase 
of SrFeAsF is also investigated. Based on the Griineisen approximation and Debye approximation for internal 
energy, the volume of SrFeAsF can be well fitted with Debye temperature of 347(5) K. The experimental atomic 
displacement parameters for different crystallographic sites in SrFeAsF are analyzed with Debye model. The 
results suggested that the expansion of FeAs layers plays an important role in determining the thermal expansion 
coefficient. 



PACS numbers: 74.70.Xa, 75.25.-j, 65.40.De 

High temperature superconductivity has recently been dis- 
covered in carriers doped iron pnictide compounds HHHHIil]. 
Till now, the highest T c attained is 57.4 K in the electron 
doped '1111' compound Can.4Nan.6FeAsF |5J, while for the 
so called '122' compound the highest T c of 38 K is reached 
in the hole doped Bao.6Ko.4Fe2As2 |4J]. It is generally believed 
that the superconductivity in iron pnictides is unlikely due to 
simple electron-phonon coupling, as demonstrated from ex- 
tensive studies of phonon dynamics 10,0]. Magnetism seems 
to play a crucial role in the appearance of superconductivity 
and AFM spin fluctuations have thus been suggested to be a 
possible pairing mechanism |8]. 

The spin density wave (SDW) transition, which is asso- 
ciated with long range antiferromagnetic (AFM) order of 
Fe moments, is shown in undoped FeAs-based compounds 
|9l iKj LLil LL2[]. Band structural calculations suggested that 
the SDW is driven either by Fermi surface nesting between 
electron and hole pockets or by the magnetic exchange cou- 
pling between local moments u2, 14, 15, 3l- Furthermore, 



the SDW transition is found to be accompanied by the tetrag- 
onal to orthorhombic (T-O) structural phase transition in the 
122 family a nd p receded by the T-O transition for the 1111 
family lf9l TTot fTlL fTsll . It is generally accepted that the crystal 
and magnetic structure in undoped FeAs compound are inti- 
mately coupled. Theoretical studies suggested that the struc- 
tural phase transition is driven by the AFM stripe magnetism 
ii. 



directly 



Therefore, it is important to investigate 



the relationship between structural properties and AFM order- 
ing of the undoped FeAs compound, especially for the 1111 
compound, in which the structural phase transition always oc- 



curred at higher temperature than SDW transition tempera- 
ture. 

As a new parent phase of iron pnictide family, the SrFeAsF 
also exhibited prominent anomaly in magnetic and electric 
property measurements 1 18 , T^, 231 • However, detailed stud- 
ies have not been carried out to investigate the magnetic and 
structural transformation in SrFeAsF. In this paper, we inves- 
tigated the structural and magnetic phase transition as well as 
thermal expansion in SrFeAsF compound by using the neu- 
tron powder diffraction (NPD) method. Both structural and 
magnetic phase transitions are clarified; the temperature de- 
pendence of atomic displacement parameters and thermal ex- 
pansion are well modeled based on the experimental results. 
The SrFeAsF polycrystalline sample was synthesized by a 
solid state reaction method as described in Ref 11811 with im- 
purity phases (CaF2 and Fe2C>3) of less than 1%. The neu- 
tron powder diffraction measurements were performed on the 
high flux powder diffractometer D20 at Institut Laue Langevin 
(Grenoble, France). A Ge (115) monochromator was used to 
produce a monochromatic neutron beam of wavelength 1.88 
A. The sample was loaded in a vanadium sample holder and 
then installed in the liquid helium cryostat that can generate 
temperature down to 2 K. The program FULLPROF 112 ill was 
used for the Rietveld refinement of the crystal and the mag- 
netic structures of the compounds. 

Similar to the CaFeAsF compound, the SrFeAsF also crys- 
tallized in orthorhombic structure with space group Cmma at 
2 K. Refinement of the neutron diffraction pattern gave lat- 
tice parameters of a = 5.6689(1) A, b = 5.6260(1) A and c 
= 8.9325(2) A at 2 K. Both lattice parameters and unit-cell 
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FIG. 1: (Color online) (a)Neutron powder diffraction pattern of Sr- 
FeAsF at 2 K. The cirlces represent the observed intensities, the solid 
line is the calculated pattern. The difference between the observed 
and calculated intensities is shown at the bottom. The vertical bars 
indicate the expected nuclear and magnetic Bragg reflection posi- 
tions, (b) The (103),k magnetic reflections obtained by subtracting 
the NPD pattern measured at 240 K from the pattern measured at 2 
K. (c) Temperature dependence of the integrated intensity of (103) M 
magnetic Bragg reflection. The solid curve is fit to the power law. 



volume V are larger than that in CaFeAsF, which is consistent 
with the fact that the atomic radius of Sr is larger than that 
of Ca. The refined neutron powder diffraction pattern of Sr- 
FeAsF at 2 K is shown in Fig. 1(a). Apart from the nuclear 
Bragg reflections, the magnetic peaks can be fitted by using 
the same AFM structural model as in CaFeAsF [22] for Fe 
magnetism and the Fe moment is deduced to be 0.58(6) \ib. It 
seems that the AFM configuration of Fe spins and the small 
magnitude of moment are a common feature for all FeAs sys- 
tems. The origin of that small iron moment in these com- 
pounds was explained theoretically as the result of the itiner- 
ant character of iron spins [13] or the nearest and next nearest 
neighbor superexchange interactions between Fe ions which 
give rise to a frustrated magnetic ground state 11411 . The tem- 
perature dependence of integrated intensity of (103)m mag- 
netic reflection was plotted in Fig. 1(c) to obtain the AFM 
phase transition temperature. By fitting the ordering parame- 
ter with power law, the Fe spins are found to order below 
- 133(3) K. The onset of the structural phase transition from 
P4/nmm to Cmma takes place at T$ = 180(2) K as indicated 
by the evolution of (220)o tetragonal reflection with temper- 
ature [Fig. 2(a)]. The T$ is in accordance with the structural 
phase transition observed by other methods |18, T^, 20ll . 
The variation of lattice parameters can be obtained through 
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FIG. 2: (Color online) (a) Evolution of (220) r reflections with the 
change of temperature. The structural phase transition is clearly re- 
vealed by the splitting of tetragonal (220) r into orthorhombic (400)o 
and (040)o reflections. (b)(c) Temperature dependence of the esti- 
mated lattice parameter a, b and c for both tetragonal and orthorhom- 
bic phases. 



the refinement of NPD patterns. As shown in Fig. (b) and 
(c), the lattice parameters a and b deviated each other when 
the temperature is lower than the structural phase transition 
temperature, while the lattice parameter c decreased gradually 
with decrease of temperature. The structural distortion will 
result in an inequivalent of nearest- and next-nearest-neighbor 
coupling (J i and J2) in FeAs layers and lift the magnetic frus- 
tration. The collinear AFM spin configuration is stabilized as 
the ground state when J 2 > JJ2 QEal- Within the Heisen- 
berg magnetic exchange model, it is established that the dif- 
ference between the structural and AFM transition tempera- 
ture is strongly correlated with the ratio between out-of -plane 
magnetic exchange coupling J z and J 2 [16,23]. To our knowl- 
edge, the difference between T$ and Tn in SrFeAsF (~ 47 K) 
is the largest among all known iron pnictides. Therefore, the 
big temperature difference between T s and Tn indicates the 
smallest JJJ% ratio in SrFeAsF. 

The detailed structural information for SrFeAsF at 2 K and 
240 K, as obtained from NPD data, are given in Table I. The 
variation of volume with temperature over measured temper- 
ature range is shown in Fig. 3. To model the experimental 
data of the unit cell volume, we followed the approach of 
the Griineisen approximation for the zero pressure equation 
of state, in which the effects of thermal expansion are con- 
sidered to be equivalent to elastic strain [24]. Thus, the tem- 
perature dependence of the volume can be described by V(T) 
= yU(T)IK() + Vq, where y is a Griineisen parameter, Kq is 
the incompressibility and Vo is the volume at T = K. By 
adopting the Debye approximation, the internal energy U(T) 
is given by: 



U(T) = 9Nk B T 
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TABLE I: Refined results of the crystal and magnetic struc- 
tures for SrFeAsF at 2 and 240 K. The atomic positions 
for space group Cmma: Sr(4g)(0,0.25,z), Fe(4£)(0.25,0,0.5), 
As(4g)(0,0.25,z), F(4a)(0.25,0,0); for P4/nmm: Sr(2c)(0.25,0.25,z), 
Fe(2£)(0.75,0.25,0.5), As(2c)(0.25,0.25,z), F(2a)(0.75,0.25,0). 
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where N is the number of atoms in the unit cell, kg is 
the Boltzmann's constant and 6q is the Debye temperature. 
Above model exhibits a good fit to volume variation as indi- 
cated by the solid line in Fig. 3. We obtained following phys- 
ical parameters: G D = 347(5) K, V = 284.89(2) X 3 and y/K Q 
= 1.57(3)xl0 -9 Pa -1 . The value for the Debye temperature as 
determined from neutron experiments is in reasonable agree- 
ment with those determined from specific heat measurements, 
D = 339(1) in Refill. Note that the Debye temperature ap- 
proximates only the acoustical phonons at low temperature. 
The thermal expansion coefficient a(T) which corresponds to 
the first derivative of U(T) divided by the unit cell volume at 
that temperature is plotted in the inset of Fig. 3(a) as the solid 
dots. The solid curve correspond to the calculated value of 
a(T) from V(T). 

As a basic crystallographic parameter, the atomic displace- 
ment parameter (ADP) can reflect the atomic thermal motion 
and provide useful information on the thermal properties of 



materials 11251 12611 . Therefore, the ADPs of different crys- 
tallographic sites in SrFeAsF was also investigated. For the 
low temperature Cmma phase, the space group symmetries 
of all special positions [Sr(4g), Fe(4b), As(4g) and F(4a)] al- 
low three independent ADP elements. Based on approximate 
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FIG. 3: (Color online) (a) Temperature dependence of the volume of 
unit cell of SrFeAsF. The solid line is fit of the Debye model as de- 
scribed in the text. The inset shows the thermal expansion constant 
of volume as determined from the fit of the Debye model in compar- 
ison to the experimental data. (b)(c) Isotropic thermal parameters for 
different atomic sites in SrFeAsF. 



classification, the isotropic ADP (U 2 ) was applied by assum- 
ing that the displacements are the same in all directions. The 
thermal parameter B, so which is obtained assuming isotropic 
thermal motions of atoms are evaluated using the relation B, - w 
= %ji 2 {U 2 ) + B sta , where B sla is the static component of the 
thermal parameter caused by the presence of a certain amount 
of static disorder in compounds. Fig. 3(b) and Fig. 3(c) 
show the temperature variation of experimental thermal pa- 
rameters for all four atomic sites of SrFeAsF. Debye model 
was adopted to describe the ADP as a function of tempera- 
ture: 
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where <$(6d/T) is given by 
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The fitted results for the thermal parameters based on above 
equations are plotted in Fig. 3(b) and Fig. 3(c) as the solid 
lines in comparison with the experimental results. The De- 
bye temperatures evaluated for different atomic species are 
Sr [238(6)K], Fe [365(6)K], As [348(6)K] and F [611(6)K]. 
The averaged mass weighted value for the Debye temperature 
is obtained to be 332(10) K. The corresponding Debye fre- 
quencies for Sr, Fe, As, F sites are 4.9(1), 7.6(1), 7.2(1) and 
12.7(1) THz respectively. The small difference between the 
Debye frequencies for Fe and As revealed the similar range of 
phonon frequencies for those two sites. The range of phonon 
frequencies obtained from our analysis is in very good agree- 
ment with our recent inelastic neutron scattering and ab initio 
phonon calculation results [27] for SrFeAsF indicating the re- 
liability of our analysis. 



ab initio electronic structure calculations for SrFeAsF |28] 
shows that it possessed essentially the same band dispersions 
in the vicinity of the Fermi surface as in other 1111 com- 
pounds, such as LaFeAsO and CaFeAsF. Furthermore, the 
electronic states near the Fermi surface are dominated by con- 
tributions from Fe and As, which indicates that the FeAs lay- 
ers are playing an important role in magnetism and supercon- 
ductivity of the compound. From our thermal parameter anal- 
ysis, we also notice that the mean value of the Debye temper- 
ature for Fe and As [355(8) K] in SrFeAsF agrees well with 
the Debye temperature deduced from the analysis of thermal 
expansion curve [Fig. 3(a), 9o = 347(5) K ]. Therefore, the 
thermal expansion of FeAs layers is supposed to be the main 
contribution to the thermal expansion of volume. 

In summary, the structural and magnetic phase transition as 
well as the thermal expansion of SrFeAsF was investigated by 
using the neutron diffraction method. The onset of tetragonal 
to orthorhombic phase transition is at 180 K while the para- 
magnetic to antiferromagnetic phase transition takes place at 
133 K. SrFeAsF exhibites the largest difference between T$ 
and Tn among all known iron pnictides. The striped AFM ar- 
rangement of Fe spins is confirmed by refinement of neutron 
data. The magnitude of Fe moment aligned along a direction 
is deduced to be 0.58(6) fig- The lattice volume of low tem- 
perature orthorhombic phase can be modeled well in frame 
of Griineisen approximation with a Debye approximation for 
the internal energy. The analysis on atomic displacement pa- 
rameters of different crystallographic sites suggested that the 
vibration of Fe and As is coupled in FeAs layers and they con- 
tribute mainly to the thermal expansion of volume in SrFeAsF 
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